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We investigate the exclusive rare B c — > D( a \vti>i and B — > D( s )£ + £~ (£ — e, fj,, r) decays within 
the standard model and the light-front quark model constrained by the variational principle for the 
QCD motivated effective Hamiltonian. The form factors f±(q 2 ) and fr(q 2 ) are obtained from the 
analytic continuation method in the q + — frame. While the form factors f+(q 2 ) and fr{q 2 ) are 
free from the zero-mode, the form factor /_ (q 2 ) is not free from the zero-mode in the q + = frame. 
We discuss the covariance(i.e. frame-independence) of our model calculation and quantify the zero- 
mode contributions to f-(q 2 ) for B c — > Dr 3 \ decays. The branching ratios and the longitudinal 
lepton polarization asymmetries are calculated with and without the long-distance contributions. 
Our numerical results for the non-resonant branching ratios for B c — > D(D S )^V£V£ and B c — s- 
D(D s )£ + £~ are in the order of 10~ 8 (1CT 7 ) and 10" 9 (10~ 8 ), respectively. The averaged values 
of the lepton polarization asymmetries obtained from the linear (harmonic oscillator) potential 
parameters are found to be -0.99 (-0.99) for B c -> D[i+ yT and -0.16 (-0.15) for B c -s> Dt + t~, 
and -0.98 (-0.98) for B c -> D s ^ + /i" and -0.14 (-0.12) for B c ->■ D s t + t~, respectively. 



I. INTRODUCTION 

The LHCb (Large Hadron Collider beauty) experiment 
dedicated to heavy flavor physics at LHC make precision 
tests of the standard model (SM) and beyond the SM ever 
more promising. Its primary goal is to look for indirect 
evidence of new physics in CP violation and rare decays 
of beauty and charm hadrons. Especially, a stringent test 
on the unitarity of Cabibbo-Kobayashi-Maskawa (CKM) 
mixing matrix in the SM will be made by this facility. 
With the upcoming chances that a numerous number 
of B c mesons will be produced at LHC, one might ex- 
plore the exclusive rare B c decays to (D, D s )£ + £~ and 
(D, D s )u£V£ (£ — e,/i,r) induced by the flavor-changing 
neutral current b — » (d, s) transitions. Since in the SM 
the rare B c decays are forbidden at tree level and occur 
at the lowest order only through one-loop diagrams [H-0| , 
they are well suited to test the SM and search for physics 
beyond the SM. In such exclusive rare decays, any reli- 
able extraction of the perturbative effects encoded in the 
Wilson coefficients of the effective Hamiltonian requires 
an accurate separation of the nonperturbative contribu- 
tions, which are encoded in the hadronic form factors. 
This part of the calculation is model dependent since it 
involves nonperturbative QCD. Therefore, a reliable es- 
timate of the hadronic form factors for the exclusive rare 
B c decays is very important to make correct predictions 
within and beyond the SM. 

There are some theoretical approaches to the calcula- 
tions of the exclusive rare B c — > (D. D s )£ + £~ and B c — s- 
(D, D s )vivt decay modes. Although we may not be able 
to list them all, we may note here the following works: 
the relativistic constituent quark model the light- 
front (LF) and constituent quar k model (CQM) [9j, and 
three point QCD sum rules |10|. The rare B c -> D s £ + £~ 
decay beyond the SM has also been studied in [ll| . Per- 
haps, one of the most well-suited formulations for the 
analysis of exclusive processes involving hadrons may 
be provided in the framework of light-front quantiza- 



tion [3. 

The purpose of this paper is to extend our our light- 
front quark model (LFQM) [ll-Ql] based on the QCD- 
motivated effective LF Hamiltonian to calculate the 
hadronic form factors, decay rates and the longitudinal 
lepton polarization asymmetries (LPAs) for the exclusive 
rare B c — > (D,D s )£ + £~ and B c — > {D 1 D s )v£V£ decays 
within the SM. The LPA, as another parity-violating ob- 
servable, is an important asymmetry [20j | and could be 
measured by the LHCb experiment. In particular, the r 
channel would be more accessible experimentally than e- 
or /i-channels since the LPAs in the SM are known to be 
proportional to the lepton mass. 

In our previous LFQM analysis [ill [l9j] , we have ana- 
lyzed the exclusive semileptonic B c — > (£>, r) c , B, B s )£i>i 
decays [l8| and the nonleptonic two-body decays of B c 
mesons such as B c — > (Z?( s ), r\ c , B( S ))(P, V) decays [l9| 
(here P and V denote pseudoscalar and vector mesons, 
respectively). Our LFQM [l3l - [l9j analysis compared 
to the other LFQM has several salient features: (i) 
We have implemented the variational principle to the 
QCD motivated effective LF Hamiltonian to enable us 
to analyze the meson mass spectra and to find opti- 
mized model parameters [HI, 0|- ip) The weak form 
factors f±(q 2 ) for the semileptonic decays between two 
pseudoscalar mesons are obtained in the Drell-Yan-West 
(q+ = q° + q 3 = 0) frame [H| (i.e., q 2 = -q^ < 0) 
and then analytically continued to the timelike region 
by changing to — q 2 in the form factor. The covari- 
ance (i.e., frame independence) of our model has been 
checked [18| by performing the LF calculation in the 
q + = frame in parallel with the manifestly covariant 
calculation using the exactly solvable covariant fermion 
field theory model in (3 + 1) dimensions. We also found 
the zero-mode 22] contribution to the form factor f-(q 2 ) 
and identified [18| | the zero-mode operator that is convo- 
luted with the initial and final state LF wave functions. 

Specifically, in the present analysis of exclusive rare 
B c -t (D,D s )£ + £~ and B c -t (D,D s )v£V t decays, three 
independent hadronic form factors, i.e. f+(q 2 ), f-(q 2 ) 
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from the vector-axial vector current, and /t(<Z 2 ) from the 
tensor current, are needed. While the two form factors 
/+ and fx can be obtained only from the valence con- 
tributions in the q + — frame without encountering the 
zero-mode complication, the form factor /_(q 2 ) receives 
the higher Fock state contribution (i.e., the zero mode 
in the q + = frame or the nonvalence contribution in 
the q + > frame) within the framework of LF quantiza- 
tion. Thus, it is necessary to include either the zero-mode 
contribution (if working in the q + — frame) or the non- 
valence contribution (if working in the q + > frame) 
to obtain the form factor /_(q 2 ). In this work, we shall 
use the form factors f+(q 2 ) and f-{q 2 ) for the exclusive 
semileptonic B c — > (D,D S ) decays obtained in [l8| and 
the form factor fr{q 2 ) obtained in [l6| for the analysis 
of B — > K£ + £~~ decay. Especially, the Lorentz covari- 
ance of our tensor form factor /t(<? 2 ) is discussed in this 
work. The present investigation further constrains the 
phenomenological parameters and extends the applicabil- 
ity of our LFQM [TU, [l4[ to the wider range of hadronic 
phenomena. 

The paper is organized as follows. In Sec. II, the 
SM operator basis, describing the b — > (d, s)£ + £~ and 
b — > (d, s)vivi transitions, is briefly presented. In Sec. 
Ill, we briefly describe the formulation of our LFQM 
and the procedure of fixing the model parameters using 
the variational principle for the QCD motivated effec- 
tive Hamiltonian. We present the LF covariant forms 
of the form factors f±(q 2 ) and /t(<? 2 ) obtained in the 
q + = frame. In Sec. IV, our numerical results, i.e. 
the form factors, decay rates, and the LPAs for the rare 
B c — s- (D, D s )£ + £~~ and B c —> (D, D s )vivt decays are 
presented. Summary and discussion of our main results 
follow in Sec. V. In the Appendix, we explicitly show 
the covariance of /t(<? 2 ) by performing the LF calcula- 
tion in parallel with the manifestly covariant one using 
the exactly solvable covariant fermion field theory model 
in (3 + 1) dimensions. 
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FIG. 1: Loop diagrams for B c — > D q (£ + £ ,viVi){q = d, s) 
transitions. 



ture constant, and are the CKM matrix elements. 
The relevant Wilson coefficients Ci can be found in 
Ref. 0. The effective Hamiltonian responsible for the 
b — > qv(i>f (q = d, s) decay processes is given by [23j j 
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where Xt = (m t /Mw) 2 and X(xt) is an Inami-Lim func- 
tion which is given by 
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The long distance (LD) contribution to the exclusive 
B c — > D q (q = d, s) decays is contained in the meson ma- 
trix elements of the bilinear quark currents appearing in 
Hlfi* and T-CIq 1 ■ In the matrix elements of the hadronic 
currents for B c — > D q transitions, the parts containing 
7 5 do not contribute. Considering Lorentz and parity in- 
variances, these matrix elements can be parametrized in 
terms of hadronic form factors as follows: 

J» EE {D q \q^b\B c ) = f + (q 2 )P» + f-(q 2 )q», (4) 



II. EFFECTIVE HAMILTONIAN 

In the SM, the exclusive rare B c — > D q (£ + £~, vtrDt) 
decays are at the quark level described by the loop 
b — > q(£ + £~, V(pi){q = d,s) transitions, and receive con- 
tributions from the Z(7)-penguin and VF-box diagrams 
as shown in Fig. [TJ 

The effective Hamiltonian responsible for the b — > 
q£ + £~(q = d,s) decay processes can be represented in 
terms of the Wilson coefficients, Cy ff , Cg S and C\q as 
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where Gf is the Fermi constant, a cm is the fine struc- 



and 
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where P = Pb c + PD q and q = P Bc - Po q is the 
four-momentum transfer to the lepton pair and 4m 2 < 
q 2 < (Mb c — MdJ 2 . We use the convention = 
(i/2)[7 M , 7"] for the antisymmetric tensor. Sometimes 
it is useful to express Eq. Q in terms of f+(q 2 ) and 
fo(q 2 ), which are related to the exchange of l - and + , 
respectively, and satisfy the following relations: 



/+(0) = /o(0), fo(q 2 ) = f+(q 2 ) + 
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With the help of the effective Hamiltonian in Eq. ([T]) 
and Eqs. (j3| and (JSJ) , the transition amplitude A4 = 



3 



(D q t+£-\H e &\B c ) for the B c -> D q t + i~ decay can be 
written as 
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The differential decay rate for the exclusive rare B c 
D q £ + £~ with nonzero lepton mass is given by [HI, 
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with s = q 2 /Ml , i = mj/M 
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m b /M Bc and 

r = M|, /M§ . The differential decay rate in Eq. ([5]) may 
be written in terms of (/+, fo, fx) instead of (/+, fx) 
as discussed in [l6j]. Note also from Eqs. ([8]) and © that 
the form factor /_((j 2 ) contributes only in the nonzero 
lepton (nig ^ 0) mass limit. Dividing Eq. (|8]) by the total 
width of the B c meson, one can obtain the differential 
branching ratio dBR(B c -s- D q £ + £~)/ds = (dT(B c -> 
D q l+£-)/T tot )/ds. 

The differential decay rate for B c — > D q vgvi can be 
easily obtained from the corresponding formula Eq. 
for B c — > D q £ + £~ by the replacement 



mi 



0, Cf 



o, a 



cff 



sin 2 9w 



C 



10 



X(a; t ) 



sin 6w 
( 10 ) 

where is the Weinberg angle. As another interesting 
observable, the LPA, is defined as 
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where h = + 1 ( — 1) denotes right (left) handed £ in the 
final state. From Eq. ([5]), one obtains for B c 

1/2 
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Because of the experimental difficulties of studying the 
polarizations of each lepton depending on s and the Wil- 
son coefficients, it would be better to eliminate the de- 
pendence of the LPA on s, by considering the averaged 



form over the entire kinematical region. The averaged 
LPA is defined by 



(Pl) = 
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III. 



FORM FACTORS IN LIGHT-FRONT 
QUARK MODEL 



The key idea in our LFQM [13J, |l4| for mesons is to 
treat the radial wave function as a trial function for 
the variational principle to the QCD-motivated effective 
Hamiltonian saturating the Fock state expansion by the 
constituent quark and antiquark. The QCD-motivated 
Hamiltonian for a description of the ground state meson 
mass spectra is given by 
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where k = (k_L,fc z ) is the three- momentum of the con- 
stituent quark, M qq is the mass of the meson, and 
l^nim) i s t ne meson wave function. We use two interac- 
tion potentials V qq ; (i) Coulomb plus harmonic oscillator 
(HO) and (ii) Coulomb plus linear confining potentials. 
The hyperfine interaction essential to distinguish pseu- 
doscalar and vector mesons is also included; viz., 



Vqq = Vq+ Vhyp = O + V c onf ~ 
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where V CO nf = br (r 2 ) for the linear (HO) potential and 
(S q ■ S q ) = 1/4 (—3/4) for the vector (pseudoscalar) 
meson. Using this Hamiltonian, we analyze the meson 
mass spectra and various wave-function-related observ- 
ables, such as decay constants, electromagnetic form fac- 
tors of mesons in a spacelike region, and the weak form 
factors for the exclusive semileptonic and rare decays of 
pseudoscalar mesons in the timelike region [l3l - [T8j . 

The momentum-space LF wave function of the ground 
state pseudoscalar mesons is given by 

*S8o(si,kij_,Ai) =n° x LX3 {x i ,-k i± )ct>(x h -k i± ), (16) 

where <fi(xi, kjj_) is the radial wave function and 7^\iA 2 * s 
the spin-orbit wave function. The model wave function in 
Eq. (|16[) is represented by the Lorentz-invariant internal 
variables, Xi = pi / P + , k,^ = p,^ — XiP± and A^, where 
pp. = (p+ 7 p- ; p ± ) = (P° + P 3 ,(A/ 2 + P 2 L )/P+,P ± ) 
is the momentum of the meson M , and pf and A^ are 
the momenta and the helicities of constituent quarks, re- 
spectively. The covariant forms of the spin-orbit wave 
function for pseudoscalar mesons is given by 
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where M = y/M§ - (m~i - m 2 ) 2 and M 2 = £ 2 = i ( k iJ_ + 
mf)/xi is the boost invariant meson mass square ob- 
tained from the free energies of the constituents in 
mesons. For the radial wave function 0, we use the Gaus- 
sian wave function: 

4 3/4 / QJ, 

<P( Xl: ] ll± ) = — ] J^±eM~k 2 /2l3 2 ), (18) 

where /3 is the variational parameter and yj dk z /dx is the 
Jacobian of the variable transformation {x,k^} k = 
(ki,fc,). 

We apply our variational principle to the QCD- 
motivatcd effective Hamiltonian first to evaluate the ex- 



pectation value of the central Hamiltonian Hq + Vq, i.e., 
(0|(i/o + Vo)|0) with a trial function (j>(xi,\s.i±) that de- 
pends on the variational parameter f3. Once the model 
parameters are fixed by minimizing the expectation value 
((j>\(Ho + Vq)\4>), then the mass eigenvalue of each meson 
is obtained as M q q = (4>\(H + V q q)\(j>). A more detailed 
procedure for determining the model parameters of light- 
and heavy- quark sectors can be found in our previous 
works [3 [II 01. 

The form factors f±(q 2 ) and /t(<? 2 ) for B c (qiq) — > 
P(<72<?) decays are obtained from the q + = frame. 
Although the form factors f±(q 2 ) and /r(? 2 ) are given 
in [lU and [Tf3], respectively, we list them here again: 
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where k'j_ = k^ + (1 — a;)qj_, A = (1 — i)ra t + im 
(i = 1,2), and q-P = M?-M$ with M\ and M 2 being the 
physical masses of the initial and final state mesons, re- 
spectively. The explicit covariances of f±(q 2 ) and fr{q 2 ) 
are proven in (l8j and in the appendix of this work, re- 
spectively. Since the form factors f±(q 2 ) and fr(q 2 ) in 
Eq. (fT9|) are defined in the spacelike {q 2 = — < 0) re- 
gion, we then analytically continue them to the timclike 
q 2 > region by changing to — q 2 in the form factors. 
We should note that our analytic method is reliable in 
the entire physical region of the exclusive rare decay since 
the first unitary branch point occurs just right after the 
zero-recoil point, q 2 = (Mi — Af 2 ) 2 . 

We also compare our analytic solutions with the double 
pole parametric form given by 



hid 2 
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1 - <7lS + (T 2 S £ 

where a\ and er 2 are the fitted parameters. 

IV. NUMERICAL RESULTS 

In our numerical calculations for the exclusive rare 
B c -> {D,D a )(v t v t ) and B c -> (D,D s )(£+£-) decays, we 



TABLE I: Model parameters {m q ,P) [GeV] for D,D S and B c 



mesons used in our 


analysis 


Q = 


= u and d. 




Model m q 


m s m c 


m b 


fiqc Psc 




Linear 0.22 


0.45 1.8 


5.2 


0.4679 0.5016 


0.8068 


HO 0.25 


0.48 1.8 


5.2 


0.4216 0.4686 


1.0350 



use two sets of model parameters (m g ,/3) for the linear 
and HO confining potenti als g iven in Table U jl3|, Q, EH ■ 
Although our predictions [18j of ground state heavy me- 
son masses are overall in good agreement with the experi- 
mental values, we use the experimental meson masses (25j 
in the computations of the decay widths to reduce possi- 
ble theoretical uncertainties. 

Note that in the numerical calculations we take 
(m c ,mb) = (1.8,5.2) GeV in all formulas except in the 
Wilson coefficient C| ff , where (m c ,mb iPO ie) = (1.4,4.8) 
GeV have been commonly used @. For the numerical 
values of the Wilson coefficients, we use the results given 
by Ref. 0: 



Ci = 
C 4 = 

eff 



-0.248, C 2 = 1.107, C 3 = 0.011, 
-0.026, C 5 = 0.007, C 6 = -0.031, 
-0.313, C 9 = 4.344, C w = -4.669, (21) 
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and other input parameters are |VtbV t *| = 0.039, 
\V tb V* d \ = 0.008, = 129, M w = 80.43 GeV, m t = 
171.3 GeV, s\n 2 9 w = 0.2233, and t Bc = 0.46 ps. The 
effective Wilson coefficient Cg taking into account both 
the short distance (SD) and the LD contributions from cc 
resonance states (J/ip, i/)', ■ ■ ■ ) has the following form Q 

Cf(s) = C 9 +Y SD (s) + Y LD (s), (22) 

where the explicit forms of Ysd(s) and Yld(s) can be 
found in For the LD contribution Yld{s), we in- 

clude two cc resonant states J/i(;(lS) and ip'(2S) and use 
r(J/V' -» = 5.26 x 10~ 6 GeV, M J/V , = 3.1 GeV, 

r J/v , = 87 x 10~ 6 GeV for J/ip(lS) and T(ip' -> £+l~) = 
2.12 x 10~ 6 GeV. \/,.. = 3.69 GeV, I> = 277 x 10~ 6 
GeV for i/)'(2S) ||. 

In Figs. and [SJ we show the q 2 dependences of the 
form factors f±{q 2 ) and fr(q 2 ) (solid line) for the B c — > 
D and B c — > D s transitions obtained from the (a) linear 
and (b) HO potential parameters, respectively. We also 
include the results (circles) obtained from the double pole 
form given by Eq. (|2"0|) . As one can see from Figs. [5] 
and [31 our analytic solutions are well approximated by 
the double pole form. 

The form factors at the zero-recoil point (i.e., q 2 = 
9max) correspond to the overlap integral of the initial 
and final state meson wave functions. The maximum- 
recoil point (i.e., q 2 = 0) corresponds to a final state 
meson recoiling with the maximum three-momentum 
\P D(s) | = (Mf c - M% )/2M Bc in the rest frame of the 
B c meson. For B c — » D transition, while the form factors 
at q 2 = obtained from the linear [HO] potential param- 
eters are /+(0) = 0.086 [0.079], /_(0) = -0.089 [-0.070], 
and / T (Q) = -0.120 [-0.108], the form factors at 
q 2 = q 2 max are /+(£«) = 1.129 [0.752], f-(q 2 max ) = 
-1.060 [-0.723], and /t^Lx) = -1-319 [-1.034]. 
As for the zero-mode contribution to the form factor 
/_( 9 2 ), i.e. f ZM -(q 2 ) = f-(q 2 ) - fl al (q 2 ) 0, we ob- 
tain the valence contribution to f-(q 2 ) as /Z al (0) = 
-0.096 [-0.081] and /I al (^ lax ) = -1-132 [-0.823] 
from the linear [HO] potential parameters. This es- 
timates about 7% ~ 15% zero-mode contribution to 
f-{q 2 ) for the B c -> D transition. For B c D s 
transition, while the form factors at q 2 = obtained 
from the linear [HO] potential parameters are /+(0) = 
0.120 [0.126], /_(0) = -0.113 [-0.099], and / T (0) = 
—0.163 [—0.168], the form factors at q 2 = g^ ax are 
/ + ( g L x ) = 0-992 [0.868], /_( 9l 2 nax ) = -0-988 [-0.719], 
and frilmax) = —1-342 [—1.157]. We also obtain the va- 
lence contribution to f-(q 2 ) as /I al (0) = -0.121 [-0.105] 
and fl al (ql mx ) = -1.114 [-0.815] from the linear [HO] 
potential parameters. This also estimates about 7% ~ 
15% zero-mode contribution to f-(q 2 ) for the B c — > D s 
transition. 

The form factors at q 2 = and the parameters Oi of 
the double pole form for the B c ^t D and B c —> D s tran- 
sitions are listed in Tables [IT] and IIII1 respectively, and 
compared with other theoretical results [8l-[l0l|. The dif- 



(a) B ^ ~> D (Linear potential) 

' LFQM 

— o Double pole fil 




q 2 [GeV 2 ] 



(b) B c -> D (HO potential) 

LFQM 

o Double pole fit 




q 2 [GeV 2 ] 



FIG. 2: The weak form factors (solid line) for B c — >• D tran- 
sitions obtained from the (a) linear and (b) HO potential pa- 
rameters. The circles stand for the results from the double 
pole fits. 



ferences of the form factors between the linear and HO 
potential model predictions for the B c — !• D are larger 
than those for the B c —> D s . Our predictions of the 
form factors are also rather smaller than other theoreti- 
cal model predictions [sj-TToj - The upcoming experimen- 
tal study planned at the Tevatron and at the LHC may 
distinguish these different model predictions. 

We show our results for the differential branching ra- 



6 



TABLE II: Results for form factors at q 2 = of B c — > D£ + £ /vtvt decay and parameters at denned in Eq. (|20[) . 

Model / + (0) cri a 2 /-(0) ai ai /t(0) ai ai 

Linear 0.086 -3.50 3.30 -0.089 -3.38 3.09 -0.120 -3.35 3.06 

HO 0.079 -3.20 2.81 -0.070 -3.28 2.94 -0.108 -3.18 2.77 

RCQM [8] 0.186 -3.48 1.62 -0.190 -2.44 1.54 -0.275 -2.40 1.49 

CQM [9] 0.123 -3.35 3.03 -0.130 -3.63 3.55 -0.186 -3.52 3.38 

SR [10] 0.22 -1.10 -2.48 -0.29 -0.63 -4.06 -0.27 -0.72 -3.24 



TABLE III: Results for form factors at q 2 = of B c — > D s £ + £ /viVi decay and parameters ai defined in Eq. (|20p . 



Model /+ (0) 


a\ 


C"2 


/-(0) 


a\ 


ai 


MO) 


a\ 


£T2 


Linear 0.120 


-3.32 


3.09 


-0.113 


-3.36 


3.14 


-0.163 


-3.28 


3.00 


HO 0.126 


-3.10 


2.73 


-0.099 


-3.12 


2.74 


-0.168 


-3.08 


2.69 


CQM [9J 0.167 


-3.40 


3.21 


-0.166 


-3.51 


3.38 


-0.247 


-3.41 


3.30 


SR [10] 0.16 


-1.55 


-2.80 


-0.18 


-0.77 


-6.71 


-0.19 


-1.43 


-3.06 



tios for B c -> (D,D s )Y,^m in Fig. HI B c -> D£+£~ 
in Fig. El and B c — > D s £ + £~ in Fig. [6j respectively. 
The solid (dashed) line represents the result obtained 
from the linear (HO) potential parameters. For the 
B c -> (D,D s )£+g- transitions in Figs. M and EJ the 
curves with (without) resonant shapes represent the re- 
sults with (without) the LD contribution to Cg . Al- 
though the form factor f~(q 2 ) does not contribute to 
the branching ratio in the massless lepton (£ = e or 
/i) decay, it is necessary for the heavy r decay pro- 
cess. As one can see from Figs. [5] and |6l the LD con- 
tributions clearly overwhelm the branching ratios near 
J/ip(lS) and tp'(2S) peaks, however, suitable £ + £~ in- 
variant mass cuts can separate the LD contribution from 
the SD one away from these peaks. This divides the spec- 
trum into two distinct regions [111 [26[ : (i) low-dilepton 
mass, 4m 2 < q 2 < Mj,^ — S, and (ii) high-dilepton mass, 
M 2 , ± S < q 2 < g^ lax , where S is to be matched to an 
experimental cut. 

Our predictions for the nonresonant branching ra- 
tios obtained from the linear and the HO potential pa- 
rameters are summarized in Table )LVI and compared 
with other theoretical predictions [§14101] within the SM. 
Since the amplitude B c — > (D, D s )£ + £~ is regular at 
q 2 = 0, the transitions B c — > (D, D s )e + e~ and B c — > 
(D , D s ) /i + fj,~ have almost the same decay rates, i.e. in- 
sensitive to the mass of the light lepton. The branching 
ratios with the LD contributions for B c — > (D, D s )£ + £~ 
(£ = /1, r) obtained from the linear (HO) potential pa- 
rameters are also presented in Table [V] for low- and high- 
dilepton mass regions of q 2 . 

In Figs. [3 and we show the LPAs for B -> D£+£~ 
and B —> D s £ + £~ (£ = fi, r) as a function of s, re- 
spectively, obtained from the linear (solid line) and HO 
(dashed line) potential parameters. The curves with 
(without) resonant shapes represent the results with 
(without) the LD contributions. In both figures, the 
LPAs become zero at the end point regions of s. However, 
we note that if mi = 0, the LPAs are not zero at the end 
points. As in the case of the B — > Kfi + yT decay where 



TABLE IV: Nonresonant branching ratios (in units of 10 ) 
for B c -¥ {D,D s )(u e ue) and B c -¥ {D,D s ){l + r) transitions 
compared with other theoretical model predictions within the 
SM. 



Mode 



Linear HO [8] [9] 



10 



B c -*DY;ViVi 1.31 0.81 3.28 2.74 (3.38 ±0.71) 

B c ^ D S J2 vm 39 37 73 92 (49 ± 12) 

B c ->lV>- 0.18 0.11 0.44 0.40 (0.31 ±0.06) 

B c ^Dsfi + ^ 5.4 5.1 9.7 13.3 (6.1 ± 1.5) 

B c ^Dt + t~ 0.08 0.04 0.11 0.12 (0.13 ±0.03) 

B c ^D 3 t + t- 1.4 1.3 2.2 3.7 (2.3 ±0.5) 



TABLE V: Branching ratios with the LD contributions for 
B c — > (D, D„)£ + £~ for low and high dilepton mass regions of 
q 2 [GeV 2 ] obtained from the linear (HO) potential parame- 
ters. 





Mode 


4mf <q 2 < 8.5 


13.8 < < q 2 ^ 


B c 


-> Dfi + ji~ 


6.56 (4.98) x 10 _1U 


5.00 (2.38) x 10~ 1U 


B c 


-¥ Dt+t~ 




6.68 (3.01) x 10" 10 


B c - 


-» D s (i + fj,~ 


2.68 (2.75) x 10~ 8 


0.86 (0.69) x 10~ 8 


B c - 


-> D 3 T + T- 




1.09 (1.07) x 10 _s 



P L ~ -1 away from the end point regions [161.1231. 124 127|. 
the LPAs away from the end point regions are close to 
— 1 for both B c — > and B c — > Z3 s /i + /i _ transi- 

tions. In fact, the Pl for the muon decay is insensitive 
to the form factors, e.g. our Pj, away from the end point 
regions is well approximated by (27j 



Pl 



CioR.eC,' 



cir 



|c 9 off l 2 



-1, 



(23) 



in the limit of C| ff — > from Eq. (jT2j) . It also shows that 
the LPA for the fj, dilepton channel is insensitive to the 
little variation of Cf ff as expected. On the other hand, 
the LPA for the r dilepton channel is somewhat sensitive 
to the form factors. 

The averaged values of Pl without the LD contri- 
butions obtained from the linear (HO) potential pa- 
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(b) B c -> D s (HO potential) 



LFQM 

o Double pole fit 




10 

q 2 [GeV 2 ] 




FIG. 3: The weak form factors (solid line) for B c — > D 3 tran- 
sitions obtained from the (a) linear and (b) HO potential pa- 
rameters. The circles stand for the results from the double 
pole fits. 



FIG. 4: Differential branching ratios for B c — s- D ^ vtvt (up- 
per panel) and B c — > D 3 ^2 V(Vt (lower panel) obtained from 
the linear (solid line) and HO (dashed line) potential param- 
eters. 



V. SUMMARY AND DISCUSSION 



rameters arc (P L (B C -> Dfj+fj,-)) = -0.99 (-0.99), 
(P L (B C -> Dt+t-) = -0.16 (-0.15), (P L (B C 
D s n + ii-)) = -0.98 (-0.98) and (P L (B C -> D s t+ T -) = 
-0.14 (-0.12), respectively. 



In this work, we investigated the exclusive rare 
scmilcptonic B c — > (D : D s )vgvg and B c — > (D, D s )£ + £~ 
(£ = e, fi, t) decays within the SM, using our LFQM con- 
strained by the variational principle for the QCD mo- 
tivated effective Hamiltonian with the linear ( or HO) 
plus Coulomb interaction. Our model parameters ob- 
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FIG. 5: Differential branching ratios for B c — > D/j, + fi~ (up- 
per panel) and B c — > Dt + t~ (lower panel) obtained from 
the linear (solid line) and HO (dashed line) potential param- 
eters. The curves with (without) resonant shapes represent 
the results with (without) the LD contributions. 



taincd from the variational principle uniquely determine 
the physical quantities related to the above processes. 
This approach can establish the broader applicability of 
our LFQM to the wider range of hadronic phenomena. 
For instance, our LFQM has been tested extensively in 
the spacelike processes [l3[ as well as in the timelike ex- 
clusive processes such as semileptonic [3, [H, EH , rare 



FIG. 6: Differential branching ratios for B c D s )i + n~ (up- 
per panel) and B c — > D s t + t~ (lower panel) obtained from 
the linear (solid line) and HO (dashed line) potential param- 
eters. 



semileptonic [16| . radiative [17] . and nonleptonic two- 
body [lj| decays of pseudoscalar and vector mesons. 

The weak form factors f±(q 2 ) and /t(<7 2 ) for the rare 
semileptonic decays between two pseudoscalar mesons 
are obtained in the q + = frame (q 2 = — < 0) 
and then analytically continued to the timelike region by 
changing q^_ to — q 2 in the form factor. The covariance 
(i.e., frame independence) of our model has been checked 



Q 




0.1 0.2 0.3 0.4 0.5 




0.1 0.2 0.3 0.4 0.5 



0.1 



0.05 



-0.05 



-0.15 



-0.2 1 — 
0.3 



B ~>Dx x 




0.1 



0.05 



-0.05 



-0.15 



-0.2 1 — 
0.3 



B — > D x T 




0.35 



0.4 



0.45 



0.5 



FIG. 7: Longitudinal lepton polarization asymmetries for 
B c — > Dfi + n~ (upper panel) and B c — 5> Dt + t~ (lower panel) 
obtained from the linear (solid line) and HO (dashed line) 
potential parameters. 



FIG. 8: Longitudinal lepton polarization asymmetries for 
B c — > D s fi + n~ (upper panel) and B c — > D s t + t~ (lower panel) 
obtained from the linear (solid line) and HO (dashed line) po- 
tential parameters. 



by performing the LF calculation in the q + = frame in 
parallel with the manifestly covariant calculation using 
the exactly solvable covariant fermion field theory model 
in (3 + l)-dimcnsions. While the form factors f+(q 2 ) and 
fr(q 2 ) are immune to the zero modes, the form factor 
f-(q 2 ) is not free from the zero mode. Our numerical 
results show that the zero-mode contribution to the form 
factor f-(q 2 ) amounts to 7% <~ 15% for both B c ->• D 



and B c — > D s decays. 

Using the solutions of the weak form factors obtained 
from the q + = frame, we calculated the branching 
ratios for B c (D,D s )v e i7 e and B c (D,D s )£+£" 
and the LPAs for B c {D,D s )£ + £- including both 
short- and long-distance contributions from the QCD 
Wilson coefficients. Our numerical results for the non- 
resonant branching ratios for B c — > D(D s)^^ 1 ^ an d 
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B c -> D(D s )£+£- are in the order of 10~ 8 (10~ 7 ) and 
10~ 9 (10 -8 ), respectively. The averaged values of the 
LPAs obtained from the linear (HO) potential parame- 
ters are found to be —0.99 (—0.99) for B c — > D/j, + /i - and 
-0.16 (-0.15) for B c -> Dt+t', and -0.98 (-0.98) for 
B c -> D s n + ir and -0.14 (-0.12) for B c -> D s t+t-, re- 
spectively. These polarization asymmetries provide valu- 
able information on the flavor changing loop effects in the 
SM. Although the q 2 dependent behaviors of our form 
factors for the B c — > D( d s ^ transitions are not much dif- 
ferent from those of other theoretical predictions (8Hl0| , 
our results for the decay rates are slightly less than those 
of Refs. This difference essentially comes from the 

different values of the form factors at the maximum re- 
coil point and may be tested by future experiments. The 
decay rates for the B c —> Di + i~ and the LPAs for the 
B c — > D s t + t~ are also quite sensitive to the choice of po- 
tential within our LFQM. From the future experimental 
data on these sensitive processes, one may obtain more 
realistic information on the potential between quark and 
antiquark in the heavy meson system. 



Following the same procedure as in [18j for the calcula- 
tion of the form factors f±{q 2 ), we obtain the manifestly 
covariant form factor fr(q 2 ) as follows: 



N(M 1 + M 2 ) 



1 pl—X 

dx I dy 



8 7 r 2 (A?-m?)(Af-mf)./ JQ 
x [xmi + ym2 + (1 — x — y)m]C, (A. 4) 

where N = gig 2 A 2 A 2 and C is given by Eq. (14) of 
Ref. 

Performing the LF calculation of Eq. (|A.1|) in the q + = 
frame, we use the plus component of the currents to 
obtain the form factor /t(<7 2 ), i-e., 



tt F (q 2 ) = (M 1 + M 2 ) 



2q 2 P+ 



(A.5) 



The LF calculation for S£ in (fO|) can be separated into 
the on-mass shell ( p~ = p~ n ) propagating part S£ n and 
the (off-mass shell) instantaneous one using the fol- 
lowing identity: 
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f> + m = (^o„ + m) + -jl + ip -Pon)- 
Then the trace term S~ in Eq. (|A.3[) is given by 



CM CM _i_ CA 1 

T ~ J Ton °Tinst> 



(A.6) 



(A.7) 



Appendix: LF covariance of tensor form factor /t(? 2 



In the solvable model, based on the covariant Bethe- 
Salpete r (BS) model of (3 + l)-dimensional fermion field 
theory [28I43C 

rent (see Eq. ([5])) is given by 



the matrix element of the tensor cur- 



4 = %S2A?A* J 



CP 



) 4 N Al NxN k N 2 N A2 ' 



(A.l) 



where g\ and g 2 are the normalization factors which can 
be fixed by requiring both charge form factors of pseu- 
doscalar mesons to be unity at zero momentum transfer, 
respectively. The denominators in Eq. (IA.1I) . are given 
by 



N k = k 



Pi 



A, 



ie, Nj =p 2 J — m 2 
-it 0' = 1,2), 



(A.2) 



where mi, m, and m 2 are the masses of the constituents 
carrying the intermediate four-momenta p\ = P\ — k, 
k, and p 2 — P2 — k, respectively. Ai and A2 play the 
role of momentum cut-offs similar to the Pauli-Villars 
regularization 28] . The trace term Sj, is given by 



where the on-mass shell propagating part Sj- has the 
same form as S£ in Eq. (|A.3|) but with p~ — p~ n . The 
instantaneous part S^ inst is given by 



an 

°Tinst 



-2g^ + Ap^[m(p 2on ■ q) + m 2 (k on ■ q)] 
+ 2g tJ ' + Ap-[m(p loD ■ q) + m^L • g)] 
+ 2g M+ A/c _ [m 2 (pion ■ q) - mi(p 2on ■ q)}, 

(A.8) 

where Ap~ = p~ —p^ on - By doing the integration over k~ 
111 Eq. (TATl) . one finds the two LF time-ordered contri- 
butions to the residue calculations corresponding to the 
two poles in k~ , the LF valence contribution defined in 
< k + < P 2 region and the nonvalence contribution 
defined in P 2 < k + < P^~ region. The nonvalence con- 
tribution in the q + > frame corresponds to the zero 
mode (if it exists) in the q + — > limit. As we have 
shown in [l8|, [28[ , the LF valence contribution comes ex- 
clusively from the on-mass shell propagating part and 
the zero mode from the instantaneous part. This implies 
that the form factor /t(<7 2 ) is free from the zero mode 
since 

Srinst = °- The LF form factor /t F (9 2 ) is then 
obtained as 



= Tr[ 75 (^i + m x )io» v q v {i> 2 + m 2 ) 75 (- # + m)] fr (« ) 
= -4{p^[m(p 2 • q) + m 2 (k ■ q)] - p 2 \m{v\ ' o) 
+mi(fc ■ q)} + fc"[mi(p 2 ■ q) - m 2 ( Pl ■ g)]}.(A.3) 



(mi — m 2 ) 



J rf 2 k±xi(x,k ± ) 



qi 



Ai 



(A.9) 
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where k'j^ = kj^ + (1 — x)q±, Ai = (1 — x)rrii + xm(i = 
1, 2), and q ■ P = M\ - Aff. The LF vertex functions xi 
and \2 are given by 



Xi(^ k ±) 



1 



x 2 (M? - Ml ) (M? -M\y 
1 

x\Ml-M'i){Ml-M'l) 



,(A.10) 



where 



M 2 



+ m 2 
1 - x 



m 



k' 2 



(A.11) 



1 — x x 
and M Ai = M 2 (mi -> A x ), Af£ = M^ 2 (m 2 -> A 2 ). We 



numerically confirm that our LF form factor f^ F (q 2 ) is 
exactly the same as the manifestly covariant form factor 
/r° v (9 2 )- This proves that fT F {q 2 ) is immune to the zero 
mode. 

Following the same procedure [l8| to obtain the form 
factors f±(q 2 ) within our LFQM, the form factor /t(<? 2 ) 
given by Eq. (fT9|) is obtained by the following relations: 



/27V 
^2N 



Xl(^> k -L 

1-x 

X2(x,k'_L 



2 (a;,k'x) 



1 



yl 2 



(A.12) 
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